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Figure 1: The armadillo man model (left) and the results of traditional MCF (top) compared to the results of our modified MCE Both sets of 
results show the computed surface after 2, 5, 10, and 25 semi-implicit time-steps. The zoom-in on the left hand show that our modified flow 
successfully avoids forming neck pinching singularities that commonly occur in the mean-curvature flow of non-convex regions. 



Abstract 

This work considers the question of whether mean-curvature flow 
can be modified to avoid the formation of singularities. We analyze 
the finite-elements discretization and demonstrate why the original 
flow can result in numerical instability due to division by zero. We 
propose a variation on the flow that removes the numerical instabil- 
ity in the discretization and show that this modification results in a 
simpler expression for both the discretized and continuous formu- 
lations. We discuss the properties of the modified flow and present 
empirical evidence that not only does it define a stable surface evo- 
lution for genus-zero surfaces, but that the evolution converges to a 
conformal parameterization of the surface onto the sphere. 

CR Categories: 1.3.5 [Computer Graphics]: Computational Ge- 
ometry and Object Modeling — Curve, surface, solid, and object 
representations 

Keywords: Surface Evolution, Mesh Fairing, Minimal Surfaces 

1 Introduction 

Mean-Curvature flow is one of the more basic flows that has been 
used to evolve surface geometry. It can be equivalently formu- 
lated as a flow that either (1) minimizes the gradient of the sur- 
face embedding or (2) minimizes surface area. As the former, it 
has played an essential role in the ar ea of mesh f airing, removing 
noise by smoothing the embedding |Taubin 1995| Desbrun et al.| 
[1999] . As the latter, it has been essential in the study of minimal 
surfaces |Chopp 1993]|Pinkall and Polthier 1993| . 

Despite its pervasiveness, the utility of mean-curvature flow has 
been restricted by the formation of singularities during the course 
of the flow. As a result, convergence proofs have been limited to a 
class of simple shapes (e.g. mean-convex surfaces) [Huisken 1984J . 
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In this work, we propose a modification of traditional mean- 
curvature flow that may address this limitation. We proceed by 
analyzing the finite-elements discretization of the flow, identify a 
potential for division-by-zero in the definition of the system ma- 
trix, and show that a minor modification to the flow removes the 
numerical instability, providing simpler expressions for both the 
discrete and continuous formulations of the flow. An analysis of 
this flow shows that the modification does not change the flow in 
spherical regions, and it slows down the evolution in cylindrical re- 
gions, avoiding the formation of undesirable neck pinches. Though 
we do not have a proof that it will always happen, we show numer- 
ous examples of the flow on highly concave, genus-zero, surfaces, 
demonstrating that not only is the flow non- singular, but it also con- 
verges to a conformal parameterization of the surface onto a sphere. 

An example of the modified mean-curvature flow can be seen in 
Figure [T] The original Armadillo-Man model is shown on the left. 
The top row shows the results of traditional mean-curvature flow 
(with vertices of collapsed regions fixed and taken out of the linear 
system, following Au et al. [2008 1) while the bottom row shows 
results of the modified flow. Although both approaches smooth the 
geometry, the concavity of the shape results in singularities when 
evolving with the traditional flow. For example, the (cylindrical) 
fingers collapse within 5 iterations, forming neck pinches before 
they can be merged into the hand. In contrast, the modified flow 
avoids these types of singularities, slowing down the inward flow 
of the extremities and allowing them to merge into the appendages 
before they have an opportunity to collapse. The limit behavior 
of the flow can be seen in Figure [To| demonstrating that the flow 
evolves to a conformal parameterization of the model onto a sphere. 

Outline 

The remainder of this paper is structured as follows. We present a 
review of both the continuous formulation of mean-curvature flow, 
and its discretization using finite-elements in Section|2] We analyze 
the numerical stability of the discretization in Section |3] and show 
how the flow can be modified to avoid potential division-by-zero. 
We evaluate our approach in Section |4j where we show the flow 



for a variety of genus-zero surfaces and discuss its properties. We 
conclude in Section p] summarizing our work. 

2 Review of Mean-Curvature Flow 

In this section, we briefly review mean-curvature flow. We start 
with the continuous formulation and then describe a semi-implicit, 
finite-elements discretization. Alt hough both are classic deriva- 
tions (see e.g. |Mantegazza 2011) for a modern treatment of the 
continuous formulation, and |Dziuk 199 0] for the finite-elements 
discretization), we repeat them here for completeness, as they are 
required for understanding the motivation for our modified flow. 

2.1 Continuous Formulation 

Informally, mean-curvature flow can be thought of as a flow that 
pushes a point on a surface towards the average position of its 
neighbors. As with image filtering (when replacing a pixel's value 
by the average value of its neighbors) this flow has the effect of 
smoothing out the geometry. 

Definition (MCF): Let M be a two dimensional manifold, let ^t • 
M ^ R^ be a smooth family of immersions, and let g?(-, ■) be the 
metric induced by the immersion at time t. We say that <^t is a 
solution to the mean-curvature flow if: 



dt 



^^t^t 



(1) 



where Aj 
the metric gt 



it is the Laplace-Beltrami operator defined with respect to 
■icg?. 



2.2 Finite Elements Discretization 

To model the flow in practice, we need to discretize the differ- 
ential equation. We transform the continuous system of equa- 
tions into a finite-dimensional system by choosing a function basis 
{5i, . ..,Bn} : M ^ R. Using this basis, we represent the map at 
time t by the coefficient vector x{t) = {xi{t), ...,XN{t)} cM? so 
that: 

Mp) = t^iit)Biip). 

i=l 

Although we cannot solve Equation [T] exactly, since the solution is 
not guaranteed to be within the span of the {5/}, we can solve the 
system in a least-squares sense using the Galerkin formulation: 

J (^•B,^diJt = J {At^fBi)d!it y\<i<N 

M M 

with djit the volume form induced by the metric gt. 

Setting D^ and V to be the mass and stiffness matrices of the em- 
bedding at time f [^ 

D\j = j{Bi ■ Bj)d^t L\j = -jgt i^tBi^VtBj) d^t 



M 



M 



a semi-implicit time discretization defines the linear system relating 
coefficients at time t + 3 io coefficients at time t\ 

{D'-dL')x{t^d)=D'x{t). 

In the above equations, V^ denotes the gradient operator with re- 
spect to the metric g^ 



^ We implicitly assume that either M is water-tight or that its boundary is 
fixed throughout the course of the flow. 



3 Modifying tlie Flow 

We begin this section by considering how numerical instabilities 
can arise with traditional mean-curvature flow and then propose a 
modified flow that resolves this problem. 

3.1 Numerical Instability 

A challenge of using mean-curvature flow becomes apparent when 
we compute the coefficients of the matrix D^ — 3L! by integrating 
with respect to the metric defined by the original embedding, <I>o, 
rather than the current embedding, <I>^ 

To this end, we consider how the geometry is stretched over the 
course of the flow, characterized by the endomorphism g~^ ' Str\ 
This operator is self-adjoint with respect to both g^ and gt. Its 
eigenvectors, vi and V2, define the principal directions of stretch 
(orthogonal with respect to both go and gt) and its eigenvalues, X^ 
and AI , define the (squares of the) magnitudes of stretch along these 
directions. 

Note that both the stretch directions, v/, and the stretch factors. A/, 
depend on the time parameter t. However, we omit it in our notation 
for simplicity. 

The Mass Matrix 

Using the chain rule, we obtain an expression for the (/,7)-th coef- 
ficients of the mass matrix as: 



D\j = jBi-Bj-^/\g^\dl^o. 



where ^/\gJgt\ gives the ratio of area elements. Since mean- 
curvature flow is area minimizing, the values of \gQgt\ tend to be 
small so the computation of Dj • is numerically stable. 

The Stiffness Matrix 

The situation gets more complicated when computing the stiff- 
ness matrix. The challenge here is due to the fact that, as area 
shrinks, the corresponding derivatives grow. As a result, since 
mean-curvature flow is area minimizing, there is potential for the 
Laplacian to blow up. 

To make this explicit, we decompose the gradient in the integrand 
of Lj • into orthogonal components along the principal directions of 
stretch, vi and V2. Then, using the fact that shrinking the domain 
of a function by a factor of A scales its derivative by the same fac- 
tor, the expression for the (/,7)-th coefficient of the stiffness matrix 
becomes: 



L'i 
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^l dvk dvk goivk^Vk) hh 
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(2) 



where gf ' is the {k, /)-th coefficient of the inverse of the 2x2 matrix 
whose (/,7)-th entry is gt(vi,Vj). 



Here we identify gt with the map from the tangent space to its dual so 
that the map g^ ^ ■ gt is a. well-defined map from the tangent space to itself 



Examining Equation |2] we observe that while the value 
^ -T^ —^ — V is stable, as it only depends on the values of the 
partial derivatives of the Bi along directions that are unit-length un- 
der the initial metric, go. the stretch ratios A1A2/A2A2 = A1/A2 and 
A1A2/A1A1 = A2/A1 might not be. In particular, as the stretching 
becomes more anisotropic (i.e. the mapping <^t becomes less and 
less conformal with respect to the metric go) one of the two ratios 
will tend to infinity and the integrand will blow up. 

3.2 Conformalizing the Metric 

Since the instability of the stiffness matrix is the result of 
anisotropic stretching, we address this problem by replacing the 
metric gt in the computation of the system coefficients by the clos- 
est metric that is conformal to go- 





MCF 


Heat Flow 


cMCF 


Catenoid 


r{t) = 1 


r{t) = 1 


r{t) = 1 


Sphere 


r(0 = \/l-4r 


r{t)=e-^' 


r(t) = ^/l-4t 


Cylinder 


r(0 = Vl-2r 


r{t)=e- 


r{t) = l-t 



Table 1: Radii of the catenoid, sphere, and infinite cylinder under 
the different flows, as a function of time t. 

4 Results and Discussion 

We begin by examining some examples of the conformalized mean- 
curvature flow and then proceed to a discussion of its properties. In 
our discussion, we assume that an initial embedding <I>o : M ^ R^ 
is given, and we take h to be the metric induced by this embedding, 

h = 2a. 
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4.1 Flowing Surfaces 



FEM Discretization 

For the finite-elements discretization, the implementation of this 
modification is trivial. Instead of requiring that the stiffness matrix 
be computed anew at each time- step (as is required by traditional 
mean-curvature flow), the modified flow simply re-uses the stiffness 
matrix from time t = 0. 

Specifically, since the conformalized metric gt has the same deter- 
minant as the old metric gt, the coefficients of the mass matrix are 
the same regardless of which of the two metrics we use. However, 
because the two eigenvalues of the conformalized metric are equal, 
Ai = A2, the coefficients of the modified stiffness matrix become: 
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making it independent of time t. 

Continuous Formulation 

The conformalization of the metric also results in a simpler con- 
tinuous formulation, allowing us to replace the PDE in Equation [T] 
with the following. 

Definition (cMCF): Let (M,/?) be a two dimensional Riemannian 
manifold (with metric h), let <I>^ : M ^ R^ be a smooth family of 
immersions, and let g?(-, ■) be the metric induced by the immersion 
at time t. We say that <^t is a solution to the conformalized mean- 
curvature flow if: 



^ = VWh\\^t 



(3) 



where A^ is the Laplace-Beltrami operator defined with respect to 
the metric h. (For a derivation, see Appendix [A]) 

Note that if the mapping <^t is conformal with respect to h, the 
Laplace-Beltrami operators A^ and A^ are related by: 



and the flows in Equations [T] and [3] are the same. 



To better understand how the conformalized flow evolves the em- 
bedding, we compare with two other flows. The first is traditional 
mean-curvature flow, which updates both the mass- and stiffness- 
matrix at each time-step. The second is the simple heat flow with 
respect to the metric h = go that keeps both matrices fixed, and has 
been proposed for efficient short-term flows in the context of mesh- 
fairing applications |Desbrun et al. 1999) . 

Analytic Flow 

We start by considering three simple examples for which an analytic 
expression of the flow can be computed. Understanding these sim- 
ple cases provides some intuition as to why the modified flow might 
be free of singularities. These examples include the (hyperbolic) 
catenoid, the (elliptic) sphere, and the (parabolic) infinite cylinder. 
For each of these geometries and all three flows, the evolved sur- 
faces can be characterized by their radius, r{t), as described in Ta- 
ble [T] (For a derivation, see Appendix [B]) 

Examining this table, we make several observations. 

First, for each shape, the derivatives of the three flows at ? = 
are equal, since all three start by flowing the embedding along the 
normal direction, with speed equal to the negative of the mean- 
curvature. 

Second, the catenoid remains fixed under all three flows since its 
mean-curvature is everywhere zero. 

Third, the heat flow remains stable for all shapes. This is expected 
since computing the flow is equivalent to repeatedly multiplying by 
the inverse of the system matrix (giving the characteristic exponent 
in the radius function) so the long term-behavior can be computed 
by projecting onto the lower eigenvectors of the system. 

Fourth, for the case of the sphere, both MCF and cMCF give the 
same result. This is because the mean-curvature flow of the embed- 
ding of the sphere is conformal, so we have gt = gt, and both flows 
define the same linear system. 

Finally, for the cylinder, cMCF slows the rate of shrinking so that 
the flow towards the cylinder's axis no longer accelerates with time. 
This is demonstrated in Figure [2] where the first few iterations of 
MCF and cMCF are shown. We can see that flowing with cMCF, 
the cylindrical center collapses more slowly, allowing the spheri- 
cal top to "catch-up" avoiding the formation of the singularity that 
appears in MCF. 

As we will see next, the same effect allows surface extremities to 
collapse into the main body without forming neck pinches and, for 




of this flow more formally by projecting the embedding function 
onto the lower-frequency eigenvectors of the Laplace-Beltrami op- 
erator.) 



u* 



Figure 2: Evolution of an embedding of a rounded cylinder using 
MCF (top) and the modified flow (bottom). Note that with cMCF 
the cylinder collapses more slowly, allowing the spherical cap to 
"catch-up", avoiding the singularity. 



genus-zero surfaces, evolves into a conformal parameterization of 
the surface onto a sphere. 

Empirical Evaluation 

We ran the modified flow on a number of genus-zero models. For 
each one, we defined the mass- and stiffness-matrices using the hat 
basis | ,Dziukl988J : 



A 



1 T-^- 1 + 1 T-^ 1 
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where N{i) are the indices of the vertices adjacent to vertex /, 7]y 
and 7^-^ are the two triangles sharing edge (/, j), and /3/- and p^: are 
the two angles opposite edge (ij). 

We performed the semi-implicit time- stepping using a direct 
CHOLMOD solver | ,Davis and Hager 199 9], running for 512 time- 
steps, term inating early if numerical instabilities were identified|j 
Following [Huisken 1984 |, we uniformly scaled the map after each 
step to obtain a surface with unit area. (This is equivalent to reduc- 
ing the time-step size at larger values of t, providing a finer-grained 
sampling of the flow when the surface evolves more quickly.) 

Our visualizations show the results of traditional mean-curvature 
flow (top), the heat flow (middle), and the modified flow (bottom). 
They show the original model on the left and the results of the flow 
for 2^, . . . , 2^ time-steps to the right, with per- vertex colors assigned 
using the normals from the original surface. 

Figure [3] shows an example of the flow for a dumbbell shape. Be- 
cause of the concavity at the center, traditional mean-curvature 
flow quickly creates a singularity (before the 16-th time-step with 
5 = 10~^) and the flow cannot proceed. In contrast, the modified 
flow slows down the collapse near the center, allowing the extremi- 
ties to collapse more quickly, evolving the embedding into a narrow 
ellipse which then flows back to a map onto the sphere. 

While the heat flow remains stable, it does not tend to evolve to- 
wards a smooth shape. (We have validated the long-term behavior 



Figures |4|8| show challenging examples of flow for highly non- 
convex shapes including brain gray-matter, the Lucy model, an oc- 
topus, a twisted dumbbell, and the se lf-intersecting su rface of an 
eigenfunction of the Dirac operator | Crane et al. 2011 1. Again, we 
find that MCF quickly creates a singularity and the flow cannot pro- 
ceed. Similarly, the heat flow does not evolve towards a smooth 
embedding. It is only using cMCF that the embeddings evolve to 
maps to a sphere. 

4.2 Discussion 

Several questions arise when considering the evolution given by the 
modified (area-normalized) mean-curvature flow. 

Does it converge? 

While this work is motivated by the goal of developing a variation 
of mean-curvature flow that is non-singular and converges when ap- 
plied to embeddings of genus-zero surfaces, we have only been able 
to provide experimental confirmation of this property and leave the 
proof (or the existence of a counter-example) as an open question. 
Figure |9] (left) empirically confirms the convergence of the flow for 
the models in Figure |3|7| giving the magnitude of the difference 
between successive maps at each iteration. 

What does it converge to? 



Looking at Figures |3|8| we observe that modified mean-curvature 
flow appears to always converge to a map onto the sphere. This is 
confirmed empirically in Figure |9] (right), which plots the variance 
of the distance of the mesh vertices from its barycenter, as a func- 
tion of the number of iterations. The plots shows that even though 
the flow may initially make the embedding less spherical, the vari- 
ance decays in the limit. 

How does it converge? 

We are also interested in characterizing the mapping from the orig- 
inal surface to the fixed point of the flow. Comparing the triangula- 
tion on the original surface to the triangulation on the limit surface 
(Figure[TO|, we see that the limit surface appears to preserve the as- 
pect ratio of the triangles, suggesting that the mapping to the limit 
surface is conformal. 

We confirm this empirically by measuring the quasi-conformal 
error, computed as the area-weighted average of the ratios of 
the largest to smallest singular values of the mapping's Jaco- 
bian [Sander et al. 2001 1. For the model in Figure [To] the av- 
erage quasi-conformal error is 1.034, which is comparable to the 
error of 1.033 for the conformal spherical parameterization of 
Springborn et al. |2008 1. 
Visualizations of the er- 
rors for both of these 
maps are shown in the 
inset on the right, high- 
lighting the fact that 
the quasi-conformal er- 
rors are similarly dis- 
tributed over the surface. 
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[Springborn et al. 2008] Modified iVlCF 
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^Numerical instability was defined by failure of CHOLMOD to produce 
a solution, due to the fact that linear system was not positive definite. 



More generally, Figure|9j(middle) shows the plots of these ratios for 
the shapes in Figures pfT] as a function of the number of iterations. 
We see that although the flow is not conformal, since the in-between 






Figure 3: Evolution of a dumbbell model with step-size 5 = 1x10 
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Figure 4: Evolution of brain gray -matter with step-size 3 = 5x 10 
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Figure 5: Evolution of the Lucy model with step-size 5 = 5x10 
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Figure 6: Evolution of an octopus model with step-size 5 = 2x 10 
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Figure 7: Evolution of a knot with step-size 5 = 5x10 ^. 



Figure 8: Evolution of the surface of an eigenfunction of the Dirac operator with step-size 5 = 5x10 ^. 
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Figure 10: The original armadillo-man model (left) and the sur- 
face obtained as the limit of the modified mean-curvature fiow 
(right). Zoom-ins on the triangulation show that the mapping ap- 
pears to preserve the aspect-ratio of the triangles, suggesting that 
the mapping is conformal. 



maps have a high quasi-conformal error, the evolution does appear 
to converge to a conformal map, for all the shapes. 

Proposition: 

If cMCF converges, than it converges to a map onto the sphere if 
and only if the limit map is conformal. 

Proof: 

(<^) If the map <^t* is conformal with respect to the metric h = go 
then initializing the flow with <I>o and evolving for time t* -\- s is 
equivalent to initializing the flow with <I>^s setting h = gt*, and 
evolving for time s. Thus, evolving the limit map under the mod- 
ified flow must also result in a uniform scaling of the limit map. 
Since uniform scaling is itself conformal, this implies that the limit 
map is uniformly scaled by traditional mean-curvature flow and, 
since the surface is compact, this implies that the surface is a sphere. 

(^) If <^t* is a map onto the sphere then we must have Ah<^t* = 
a^\h-^gt^\<^t* for some rescaling constant a. In particular, this 
implies that the heat-flow d^t/dt = Ah<^t will evolve <!>?* along 
directions normal to the surface <I>^*(M). Thus the function <^t^, 
considered as a map from M onto the sphere, is harmonic with re- 
spect to the metric h and therefore, by Corollary 1 of | ,Eells and| 
[Wood 1976) , conformal. 

What drives it? 

In this work, the flow was derived by modifying traditional mean- 
curvature flow. However, one can also interpret the flow as a gradi- 
ent descent on a non-negative energy. In particular, the flow is a de- 
scent on the Dirichlet energy of <^t and can be expressed as the sum 
of a (modified) area energy that drives traditional mean-curvature 
flow and a conformal energy. (For more details, see Appendix [C]) 
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Figure 9: Convergence (left), conformality (middle), and sphericity (right) of the modified me an- curvature flow for the models in Figures 3p 
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Figure 11: Evolution of a cylinder with boundary under the modi- 
fled mean-curvature flow. 
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Figure 12: Evolution of a surface with planar boundary under the 
modifled mean-curvature flow (left four images) and a bird's-eye- 
view oftheflnal mapping (right). 



Other geometries 

While the previous discussion considers embeddings of water-tight, 
genus-zero surfaces, it is also interesting to consider the behavior of 
the flow on other geometries. To this end, we have applied our flow 
both to embeddings of surfaces with boundaries, and to embeddings 
of surfaces with higher genus. 

Figures [TT] and [12] show the results of our flow for two surfaces 
with boundaries for which traditional mean-curvature flow gener- 
ates neck pinches. Though the flow converges, evaluation of the 
quasi-conformal error shows that the limit map is not conformal. 
Analyzing the open cylinder, it becomes apparent that stability and 
conformality cannot be satisfied simultaneously. If the mapping 
were conformal, then the embedding would be fixed under tradi- 
tional mean-curvature flow, implying that the corresponding surface 
is minimal. However, since the radius of the cylinder is one and the 
length of the axis is four, there is no catenoid passing through the 
two boundaries, and the only minimal surface is the one comprised 
of two disconnected disks. Thus, the mapping would only be con- 
formal if the flow were to disconnect the surface, which would be 
impossible without passing through a singularity. 

Figure [13] shows the results of our flow on two models that are not 
simply connected. Though neither flow collapses, they also do not 
appear to converge. While for higher-genus models the limit map 
cannot be conformal (as there are no embeddings of compact, non- 
genus-zero surfaces that are uniformly scaled by mean-curvature 
flow) it is possible that a limit map exists. However, even if it does, 
it is not clear that the limit is a 2-manifold. (For example, the em- 



bedding of the star appears to converge to a map onto the circle.) 

Relationsliip to 1D Flows 

While traditional mean-curvature flow of embeddings of 2D sur- 
faces in 3D can form singularities, this is not the case for em- 
beddings of ID curves in the plane. In the case of curves, the 
(uni formly rescaled ) flow always converges to a map onto the cir- 
cle [Grayson 1987| . This agrees with the empirical behavior of 
our modified mean-curvature flow in that the deformation of the ID 
curve is always conformal and the definitions of MCF and cMCF 
agree. 

Note that, as in the 2D case, (locally) scaling the map by a scales 
the Laplace-Beltrami operator by 1/a^. However, since the ID in- 
tegrals only scale by a, the two scaling terms do not cancel out and 
the discretized Laplace-Beltrami operator does not stay constant. 

5 Conclusion 

In this work, we have considered the problem of singularities that 
arise in mean-curvature flow when evolving non-convex surfaces. 
Analyzing the finite-elements discretization that commonly arises 
in geometry processing, we have associated a potential cause for 
the formation of singularities with the non-conformality of the flow. 
We have proposed a modification of the flow that simplifies both 
the discrete and continuous formulations of the flow. Although we 
do not have a proof, the work presents empirical evidence that the 
flow stably evolves genus-zero surfaces, converging to a conformal 
parameterization of the surface to a sphere. 
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A Continuous Formulation 

Although our work has focused on the finite-elements formulation 
of MCF, the modified flow can also be formulated in a continuous 
framework. 

Claim: cMCF is driven by the PDF: 



dt 



VWh\^h^t- 



Proof: To show this, we choose a test function B : M 
consider the Galerkin formulation of the PDF: 



and 



/(t-i* 



\gt'h\Ah^fByi^t. 



M 



Note that, because the metric changes over the course of the evolu- 
tion, the integrands are expressed with respect to the measure djit 
not djiii. 



Applying the change-of-coordinates formula, followed by the Di- 
vegence Theorem, we get the weak formulation: 

l(VW^\^-B]d^h = l{Ah^fB)d^h 



M 



M 



= - jh{yH^tyhB)diiH 

M 



which gives rise to the discretization with varying mass-matrix (Ihs) 
but fixed Laplace-Beltrami operator (rhs), presented in Section [3] 

B Analytic Flow Solutions 

To better understand the behavior of the flows, we consider three 
simple surfaces: The (hyperbolic) catenoid, the (elliptic) sphere, 
and the (parabolic) infinite cylinder. For each, we analyze the evo- 
lution of the embedding of the surface under the actions MCF, heat 
flow, and cMCF. 

In our analysis, we use the fact that the Laplacian of the embedding 
function is the mean-curvature weighted normal, A<I> = —HN. 

Catenoid 

Since the catenoid has mean-curvature zero, H = 0, the Laplacian 
of the embedding is zero and the surface does not evolve under any 
of the flows. 

Sphere 

Due to the rotational symmetry of the sphere, we know that its evo- 
lution under all three flows will have the form '^tip) = r{t) ■N{p) 
for some radius function r{t) and normal vector N{p) = p, so that 
d^t/dt = /{t) ■N{p). Without loss of generality, we take the initial 
radius to be one. 

Traditional MCF Using the fact that the mean-curvature of a 
sphere of radius r is 2/r, we have: 



dt 



^t^t 



IN 
r 



r(t) = Vl-4t. 



Heat Flow For the heat flow, the surface is evolved by always 
using the Laplacian at time ^ = 0: 






-- Ao<^t = r ■ Ao<I>o 



-2r-N 



r{t) = -2r{t) 
r{t) = e-'^. 



Modified l\/ICF Using Equation[3] the evolution of the surface un- 
der the modified flow can be described by scahng the heat flow by 
the reciprocal of the area change: 



d^t 


-Ir-N 


IN 
r 


=> r'{l 


^ 2 


dt 


,2 - 


)- r{t) 




) = x/l-4f. 



(Infinite) Cylinder 

Due to the translational and rotational symmetries of the cylinder, 
we know that its embedding will evolve as a constant offset from the 
original cylinder along the normal, ^t{p) = ^o{p) —N{p) + r{t) • 
N{p) for some radius function r{t), so that d^t/dt = r'{t) ■N{p). 
Without loss of generality, we take the initial radius to be one. 

Traditional MCF Since the mean-curvature of a cylinder of radius 
r is 1/r, we have: 






\^t 



N 
r 



r'{t) : 



r{t) = Vl-2t. 



Heat Flow For the heat flow, we use the fact that, on a unit-radius 
cylinder, the Laplacian of the normal is equal to the Laplacian of 
the embedding, AqA^ = Ao<I>o = -N: 



dt 



Ao(^o-N+r-N 



-r-N 



r'{t) : 
r{t)- 



-r{t) 



l\/lodified l\/ICF As above, we use Equation[3]to get: 

^ = ^^ = -iv ^ /w = -i 

dt r ^ ^ 

=^ r{t) = l-t. 



C Energy of the Flow 

We show that cMCF can be formulated as the gradient flow of an 
energy that is the sum of a smoothness term, adapted from the en- 
ergy defining traditional MCF, and a conformal energy. As we will 
see, this modified energy is just the Dirichlet energy of ^t with 
respect to the metric h. 

Energies 

Using the Euler-Lagrange formulation, traditional MCF can be de- 
fined as the gradient flow of the area functional: 



£A(^f) = Area 



M 



\h-'gt\ 



dfJ^h. 



We define a conformal energy by measuring the extent to which 
h'^gt differs from a scalar multiple of the identity. Specifically, set- 
ting Tr{h'^gt/2) ■ id. to be the scalar multiple of the identity with the 
same trace as h'^gt, we set: 



Eci^t 



I 



■-■ft-Tr(fr'ft/2)-id.||| 



|A-'ftl 



diit 



1 /Tr2(/,-'ft)-4|/!-'ft| 



:/ 



vw^\ 



djlh 



where the division by the determinant of h'^gt makes the integrand 
invariant to uniform scaling of the map O^. 

Modifying the Energies 

To define the energy driving our modified MCF, we simplify the 
energies by replacing the geometric mean of the eigenvalues in the 
denominator, \J\h-^gt\, with the arithmetic mean, Tx{h'^gt)/2, and 
sum the (modified) energies |j 

Replacing the denominators, we get: 



EAi^t 



Eci^t 



M 

1 rTr^(h-'g,)-4\h-' 



gt\ 



2 J rr{h-'g,) 

M 



dllh- 



And, taking the sum of the energies, we get: 

M 

That is, we replace the area functional defining traditional MCF 
with the Dirichlet energy of the map <^t with respect to the metric 
h. 

Linearizing the energy by considering ^ : M ^ R^ gives: 

£(cD, + £^)=^(cD,)-e|(A,cD,).WM/. + 0(£2) 

M 

= ^(cD,)-(vW^A,cD,,vp\ +0(^2), 



where (•, •)^ is the inner-product defined on the space of functions 
on M by the metric gt : 



{f,g),=l{f-g)dH, 



Thus, the gradient of the energy E, defined with respect to the inner- 
product (•, •)^, is -^\g-/h\Ah^f 

Relationship to Heat Flow 



Note that the heat flow: 






^h^t 



can also be defined as a gradient flow on the Dirichlet energy of 
O? with respect to the metric h. However, as demonstrated in Sec- 
tion [4] heat flow and cMCF evolve the surfaces in different ways. 
This is because defining the energy gradient requires choosing an 
inner-product on the space of functions on M. For heat flow, the 
inner-product is defined by the initial metric h giving (•, •)/i while 
for cMCF it is defined by the metric gt giving (•,)?• 

In particular, though both (un-normalized) flows evolve towards the 
same critical point, (the constant map taking all points in M to a 
single point in 3D), they converge to this function in different ways. 
Our modified MCF appears to converge to the constant function "as 
a sphere" (| ,Huisken 19 84]) while the heat flow does not. 

^ Note that both denominators scale quadratically with ^t, we have 
Tr{h'^gt/2) < ^\h'^gt\, and the two are equal if and only if ^t is conformal 
with respect to the metric h. 



